Evolution equations for pulse propagation in nonlinear media 



Debabrata Pal, Amitava Choudhuri and B Talukdar 
'^Department of Physics, Visva-Bharati University, Santiniketan 731235, India 

e-mail : hinoyl23@hsnl.in 



Abstract 

We show that the complex modified KdV (cmKdV) equation and 
generalized nonlinear Schrodinger (GNLS) equation belong to the Ablowitz, 
Kaup, Newell and Segur or so-called AKNS hierarchy. Both equations 
do not follow from the action principle and are nonintegrable. By in- 
troducing some auxiliary fields we obtain the variational principle for 
them and study their canonical structures. We make use of a coupled 
amplitude-phase method to solve the equations analytically and de- 
rive conditions under which they can support bright and dark solitary 
wave solutions. 



PACS numbers : 47.20.Ky, 42.81.Dp, 02.30.Jr 



1. Introduction 

In a pioneering work Zakharov and Shabat [1] solved the nonlinear Schrodinger 
(NLS) equation by the use of inverse spectral method. Subsequently, Ablowitz, 
Kaup, Newell and Segur (AKNS) [2J sought a generalization of the method 
and introduced the so-called AKNS hierarchy which involves a family of in- 
tegrable equations associated with the Zakharov-Shabat eigenvalue problem. 
It is now fairly well known that many physically important integrable partial 
differential equations belong to the AKNS hierarchy. For example, besides 
the NLS equations, this hierarchy yields the KdV-, mKdV-, Sine-Gordon-, 
Harry-Dym equations as well as the constraint KP hierarchy [3]. There are 
some relatively recent attempts to disclose the possible connections between 
the two-component Camassa-Holm equation and AKNS hierarchy extended 
by a negative flow [H [5]. Thus realization of physically important nonlin- 
ear evolution equations in the frame of AKNS model is still an interesting 
curiosity. 
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In the present work we shall provide, in Sec. 2, a derivation for two non- 
intgrable partial differential equations by using the integrable AKNS model 
for which there exists a well defined spectral problem. The equations of our 
interest are the complex modified KdV- (cmKdVI) [6] and generalized third- 
order nonlinear Schrodinger (GNLS) equations f?]. The cmKdVI equation 
describes the nonlinear steady-state propagation of lower-hybrid waves in a 
uniform plasma and more interestingly, its solutions reveal a close connec- 
tion between classical soliton and envelope solitons. On the other hand, the 
GNLS equation is used to model propagation of ultra short pulses in optical 
fibers [8] and is traditionally obtained from Maxwell's equations with special 
attention to nonlinear susceptibilities of the associated optical medium [9J. 

An awkward analytical constraint for both cmKdVI- and GNLS equations 
is that these do not follow from the action principle to have a Lagrangian 
representation which plays a role in many applicative contexts [10]. For an 
arbitrary differential equation one can always find the variational principle 
by taking recourse to the use of auxiliary field variables [11]. We follow 
this route, in Sec. 3, to construct expressions for Lagrangian densities of 
the cmKdV and GNLS equations and use them to study their canonical 
structures. In Sec. 4 we make use of a coupled amplitude-phase formulation 
[12] to solve these equations analytically and derive conditions under which 
cmKdVI- and GNLS equations can support bright and dark solitary wave 
solutions. Finally, in Sec. 5, we make some concluding remarks with a view 
to summarize our outlook on the present work. 



2. Derivation of cmKdVI- and GNLS equations using the AKNS 
hierarchy 

The AKNS hierarchy with the Zakharov-Shabat eigenvalue problem is 
given by [2] 

'(:;)-(2m(:)=o. (1) 

Here u = u{x, t) and v = v{x, t) represent two (1 + 1) dimensional complex 
fields. The operator uj{2Lzs) gives the dispersion relation of the linearization 
equation in the u component with L^g, an integrodifferential operator written 
as [13] 

~ 2. I -2v dyv -1 + 2v dyu ) ' 

The subscript zs of L merely indicates that the AKNS hierarchy is the family 
of integrable equations belonging to the Zakharov-Shabat eigenvalue prob- 
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lem. The NLS hierarchy can be obtained by choosing cc;(2L^5) = {2Lzs)"' , n = 
2, 4, 6 ... and v — —u*. In particular, for n — 2, (1) gives the NLS equation. 
For n = 4, we get the forth-order equation in the NLS hierarchy and so on. 
The choice for the dispersion relation amounts to demanding that (jj{2Lzs) is 
an entire function of the argument. 

To derive the nonintegrable equations of our interest we begin by intro- 
ducing 

oo{2L,,) = {-2L,sT (3) 

where n — 1, 2, 3.... Note that for even values of n the relation in (3) 
coincides with the dispersion relation used for the NLS hierarchy. For n — ?> 

.i2L^^)^-^(l^^ (4) 

with 

^3 rx Q px px 

On = - '2u2x J dyv - Gu^v - ^liv— + 4:u'^v J dyv + 2ux J dyvy 

r 

—2u I dyv2y , (5a) 

J — oo 

/X l-X px l-X 

dyu + 2uUx — 4:U^v / dyu + 2ux / dyVy + 2u dyu2y , 
■oo J —oo J —oo J —oo 

(56) 

/*X f'X 

0-21 = -2«a; / dyVy - 2vVx -2v dyv2y + 4«t)^ / dyv - 2v2x / dyv 



(5c) 

and 

^3 rx Q rx rx 

022 = -TTT + 2^20: / dyu + Qv^u + Auv Auv"^ / dyu-2vx / dyUy 

OX^ J-oo ox 



-\-2v I dyu2y ■ {5d) 

J —oo 

From (1), (3), (4) and (5) we get cmKdVI equation 

ut + 2{\u\'^u)^ + = (6) 
for v — —u*. Similarly, for n = 5 we get the fifth-order cmKdVI equation 

Uf — {S\u\'^U2x + 4:\u\'^U + 2u^U2x + (iUxU*)x ~ ^0u*U2x + i4:UxU2xU* — V2uUxU2x — 

lQv?U*^Ux — U5x = . (7) 
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The family of higher-order equations obtained in this way does not form a 
hierarchy in the same sense as used in the case of integrable equations. 

We have just seen that for n = 3 the dispersion relation in (3) generates 
cmKdV I equation. Let us now define 



n=2 



Written explicitly (8) becomes 



^(2L.,) = -( ?i2 + ^^i2 ) (9) 



with 



g2 rx 



/X PX 
dyv — 2u dyVy , (10a) 
-00 J —00 

/X PX 
dyu — 2u dyuy , (106) 
-00 J —00 

fX fX 

621 = —2vx I dyv — 2v dyVy (10c) 



and 



|2 rx j-x 

+ Auv + 2vx / dyu ~ 2v dyUy . (lOd) 

J ~oo J —00 



J22 



From (1), (9) and (10) we get GNLS equation 

iut + 2i\u\ Ux + 2^(1^1 )xU + U2x + 2\u\ u + iu2,x = (11) 

for V = —u*. 

3. Variational principle and canonical formulation 

The first step towards a canonical formulation of any system is to assure 
the existence of a Lagrangian. Both cmKdVI and GNLS equations in (6) and 
(11) invalidate the Helmholtz condition [Tl] such that these equations are not 
Euler-Lagrange expressions. One standard method for finding variational 
principle [11] for such equations is to introduce a set of auxiliary variables 
v{v^, ,v^} and consider the Lagrangian density C 

C{u, v) = vA[u] (12) 
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for the problem where A[m] = defines an arbitrary system of differential 

equations with A = A(Ai, ,Aj). For the cmKdVI equation in (6) we 

consider an additional complex field v — v{x, t) and introduce the Lagrangian 
density 

£, = V [ut + 4:UU*Ux + '^u^ul. + -usa:) — V* {ul + Auu*ul + 2u*'^Ux + ul^) (13) 
to write the action principle as 

5 I £dxdt = . (14) 



Clearly, the Euler-Lagrange equations for v and v* give the cmKdV I and 
its complex conjugate equations. The Euler-Lagrange equations for u and u* 
yield the equations for v and v*. The set of four equations thus obtained can 
be written in the matrix form 



/ 


u 






( 


-dl - 4:\u\'^dx 












u* 








-dl - A\u\'^dx 












V 














-dl - 4\u\'^dx 


2u*'^dx 




V 


V* 


) 


t 


{ 








2u'^dx 


-dl - 4\u\'^dx 


J 



/ u \ 

u* 

V 

V V* ) 

The canonical momentum densities are 



TT = — — = V and TT 
out 



(15) 



dC 

dul 



— —V 



(16) 



corresponding to the Lagrangian density in (13) which via Legendre trans- 
formation leads to the Hamiltonian density 



H = V* (4uw*M* + 2u*^Ux + ul^) — V (4:UU*Ux + 2w^m* + u^x) 
with the Hamiltonian written as 



H 



Hdx 



(17) 



(18) 



The equations in (15) can be written in the Hamiltonian form 

C - ^ = Mx) , H{y)} 



(19) 
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where ^ and rj stand for appropriate field variables. It is easy to verify that 
(19) is endowed with Poisson structures 

{u{x) , u{y)} = 6{x - y) , {u{x) , u*{y)} = , {u{x) , v{y)} = , 

{uix) , v*iy)} = , (20a) 
{u*{x) , u{y)} = , {u*{x) , u*{y)} = 5{x - y) , {u{x) , v{y)} = , 

{u{x),v*{y)} = , (206) 
{v{x) , u{y)} = , {v{x) , u*{y)} = , {v{x) , v{y)} = 6{x - y) , 

{v{x) , v*{y)] - (20c) 

and 

{v*{x) , uiy)} = , {v*ix) , u*iy)} = , {v*{x) , viy)} = , 

{v*{x),v*{y)} = d{x-y) . (20d) 

The Lagrangian and Hamiltonian densities for the GNLS equation in (11) 
can be written as those in (13) and (17). The Lagrangian system of equations 
is given by 

/ u\ 



1 


u 
















u 








V 






\ 


V* 







C21 C22 C23 C24 
C3I C32 C33 C34 



V 

\v* J 



(21) 



with 



and 



Cll — C44 


= -idl - dl - 4:i\u\^d^ , 


(22o) 


C12 = 




(226) 




C21 = -c*2 


(22c) 


C22 = C33 


= -io^^ + dl- 4i\u\^d^ , 


{22d) 


Cl3 = 


cu = C23 = C24 = , 


(22e) 




C31 = 2u*v , 


(22/) 




C32 = 2uv , 


(22^) 


C34 


= 2u*^ - 2iu*^d^ , 


{22h) 




* 

C4I — "^32 , 


(22z) 




C42 = -C31 


(22j) 




C43 = -C34 


{22k) 



6 



The set of equations in (21) can be written in the Hamihonian form as given 
in (19) with Poisson structures similar to those in (20) 



4. Solitary wave solutions 

To obtain the sohtary wave solutions of (6) and (11) we take recourse to 
the use of coupled amplitude-phase formulation [12] and write u{x,t) in the 
form 

u{x, t) = P{x)&^p[i{kx — ujt)] , X = X + , (23) 

with X the travelling cordinate containing the group velocity j3 of the wave 
packet. Here the function P is real. From (6) and (23) we obtain 

PP^ + 6P^Px + P^^^ - 3k^P^ + i{2kP^ -uP + 3kP^^ - k^P) = . (24) 

Equating the real and imaginary parts of (24) seperately to zero we have 

Pxxx + QP'Px + i^-^k')P^ = (25) 

and 

Pxx + Ip' - ^fp = . (26) 
The third order equation in (25) can be integrated to write 

P^x + 2P^ + 3k'^)P = . (27) 

Both equations (26) and (27) can be solved analytically. The solution of (26) 
will depend explicitly on uj and k while the solution of (27) will have similar 
dependence on f3 and k. In particular, P(x)'s will be given by 

Pix) = i^k' - P)hech[{3e - Pf^x] (26') 

or 

Pix) = (^y'^^'h[{^f)h] . (270 

Here equations (26') and (27') refer to solutions of (26) and (27) when the 
first integrals of these equations are taken as zero. Clearly, the compatibility 
condition of these solutions implies that (26') must satisfy (27) and (27') must 
satisfy (26). This viewpoint yields the necessary and sufficient condition to 
get a relation among P, u and k and we have 

^ = 2^-"^ . (28) 
k 



Since the first integral of (26) or (27) gives tlie energy E of tlie wave, we infer 
from (26') and (27') that the zero energy solution of the cmKdV equation 
represents a bright sohtory wave solution provided the velocity /? satisfies the 
constraint in (28). 



The choice 



and — ^^^ir- for the first integrals of (26) and (27) 



leads to the solutions 



P{X) = I ~ I tana 



and 



Fix)- 



CO 



2k 



tank 



u + k'' 
2k 



3\ 2 



{x + (5t) 



The compatibility condition of (26") and (27") gives the relation 



P^Ak' + - 

rv 



(26") 



(27") 



(29) 



It is of interest to note that for 13 as given in (29), the expressions for the 
first integral for (26) and (27) become identically equal. Thus we infer that 
for negative values of E the cmKdVI equation supports dark solitory wave 
solution. 

We have carried out a similar analysis for the GNLS equation in (11) and 
found that the zero energy solutions are given by 



e-k^- 

1-iO 



sech 



k^-e- 
1-cu 



{x + I3t) 



or 



P{x) = (3fc2 -2k-(3y sech {3k^ -2k-(3y{x + (3t) 
subject to the consistancy condition 

k\k-l) 



P = k{3k - 2) + 



l-cu 



+ 



CO 



(30) 



(31) 



(32) 



It is evident from (30) or (31), as with the cmKdVI equation, the zero en- 
ergy solution of the GNLS equation also represents a bright solitory wave. 
Expectedly, the negative energy solution 



Pix) 



2{l-k) 



tanh 



uj + k^ - k"" 
2(1 -3/c) 



2\ 2 



(x + (3t) 



(33) 
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or 

.2X2 



, ,'2k- p -SPY 
P{X) = ^ tank 



2k- f3-3t, , 

' + /9t) 



(34) 



with 

/5 = ^ - 2k{\ - k) (35) 



corresponds to a dark solitory wave. 



5. Conclusions 

We began by noting that a large number of physicaUy important partial 
differential equation belongs to the AKNS hierarchy. This motivated us to 
examine if the nonintegrable cmKdVI and GNLS equations could be embed- 
ded into the same hierarchy. We provided derivation of these equations in 
the frame of the AKNS model. The equations are nonLagrangian even in the 
potential representation. But we found that one can obtain the variational 
principle for them by taking recourse to the use of some auxihary fields. This 
provided a natural basis to study their canonical structures. We solved these 
equations by using a coupled amplitude-phase method and explicitly demon- 
strated the zero-energy solutions correspond to bright solitary waves while 
the negative-energy solutions represent dark solitary waves. 
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